Origin of Universal Optical Conductivity and Optical Stacking Sequence Identification 

in Multilayer Graphene 
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We show that the origin of the universal optical conductivity in a normal iV-layer graphene 
multilayer is an emergent chiral symmetry which guarantees that cr(u>) = Na un i in both low and 
high frequency limits. [a un i = (tt/2) e 2 /h]. We use this physics to relate intermediate frequency 
conductivity trends to qualitative characteristics of the multilayer stacking sequence. 
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Introduction. — Graphene is an atomically two- 
dimensional material which can be viewed either 
as a single sheet of graphite or as a large unrolled 
nanotube. Experimenters [H, d, H, 0, Q have recently 
made progress in preparing and measuring the electronic 
properties of single and multilayer graphene sheets. One 
particularly intriguing property of neutral single-layer 
graphene sheets is its interband optical conductivity 
which is expected @, 0, Q to be approximately constant 
over a broad range of frequencies with a value close to 



7r e 
27T 



(1) 



dependent only on fundamental constants of nature. 
Experiments 0, [l(| EH E3] have demonstrated that cor- 
rections to the constant universal conductivity which 
might be expected to follow from electron-electron 
interactions [Tj, EH EH or refinements of the Dirac- 
equation band-structure model[16| for example, are 
small. Recently Gaskell et al. [T3] found that for frequen- 
cies in the optical range the conductivity per layer in 
multilayer graphene sheets is also surprisingly close to 
o~ U ni- Separately Kuzmenko et aZ.[l8| demonstrated ex- 
perimentally that in bulk graphite the optical conductiv- 
ity per layer has a smoother frequency dependence and 
is even more uniformly close to cr uni than in thin mul- 
tilayers, and explained the weak frequency dependence 
they found in terms of the crossover to three-dimensions. 
In this Letter we identify the emergent chiral symmetry 
of multilayers [19|, |2fJ as a key element of the physics re- 
sponsible for the ubiquity of <7 U ni in multilayer graphene 
systems. 

A single graphene sheet consists of a honeycomb lattice 
of carbon atoms. Graphene 7r-orbitals can be viewed^ as 
possessing a "which sublattice" pseudospin degree of free- 
dom. The envelope functions of their electron waves are 
described by a two-dimensional massless Dirac equation 
which possesses pseudospin chiral symmetry and leads to 
eigenspinors in which the phase difference between sub- 
lattices 4> — ±J(f>k where the pseudospin chirality J = 1 
and cf>k = tan -1 (k y / k x ) . The pseudospin chirality is de- 
fined by this equation as rate at which pseudospin orien- 
tation varies with momentum orientation. 



When sheets are stacked to form a multilayer system 
there is an energetic preference for an arrangement in 
which each layer is rotated by 60° with respect to one 
of the two sublattices of its neighbors. This prescription 
generates three distinct planar projections of the honey- 
comb lattice (A, B, and C) and therefore 2 N ~ 2 distinct 
./V-layer sequences. We refer to multilayers in this class 
as normal. Repeated AB (Bernal) stacking and repeated 
ABC (orthorhombic) stacking should be viewed as ex- 
treme cases, as we explain in more detail below. The 
emergent chiral symmetry we discuss below applies for 
any normal multilayer. AA (hexagonal) stacking (plac- 
ing a layer directly on top of another) is energetically 
costly [2 1J, does not yield chiral symmetry, and has clear 
optical signatures. Our discussion of optical properties 
for normal multilayer graphene systems starts from a 
model, referred to below as the ideal model, in which 
only the dominant interlayer nearest-neighbor hopping 
processes couple individual-sheet Dirac-equation waves. 

The optical conductivity of an TV-layer system is ex- 
pected to approach Na un i for frequencies that exceed 
the interlayer-coupling scale but are smaller than the tt- 
bandwidth scale, since the layers then contribute inde- 
pendently and the Dirac model still applies. The ideal 
model of normal graphene multilayers has a surprising 
property which we have explained previously [l9l |20| . In 
the low-energy limit its spectrum separates asymptoti- 
cally into Nrj < N decoupled pseudospin doublets, each 
of which has chiral symmetry with a chirality J which 
can in general be larger than 1. We demonstrate below 
that the conductivity of a pseudospin doublet with chi- 
rality J is Ja un i. It then follows from the chirality sum 
rule 111], 



]T J n = N, 
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that the conductivity of the ideal model unexpectedly 
also approaches Na un i in the uj — * limit. Since the 
asymptotic pseudospins which emerge at low-energies are 
in general spread across a number of different layers and 
are qualitatively dependent on the way in which the lay- 
ers are stacked, the low-frequency limit of the interband 



2 



conductivity does not result from independent single- 
layer contributions but has a completely different origin. 
We show below that corrections to the ideal model are 
small and therefore use it to analyze correlations between 
stacking and deviations from No~ un i at intermediate fre- 
quencies. 

Chiral doublet conductivity. — The Kubo formula for the 
real part of the optical conductivity, o\r(w) = Re[cr xa; (a>)], 
of an M-band two-dimensional electron-gas system is 
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where e n j. and |n, k) are eigenvalues and eigenvectors of 
the M x M Hamiltonian matrix Tt, f n ,k is a Fermi occu- 
pation factor and v a = dH/hdk a is the velocity operator. 
The M = 2 Hamiltonian matrix of a doublet with chiral- 
ity J is 



T-tj = 7i 
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where 71 is the nearest-neighbor interlayer hopping, v^. = 
Hvke z ^ h /71 and v is the effective in-plane Fermi veloc- 
ity (for example, v — -^^p for J = 1 monolayer and 
J = 2 bilayer graphene where 70 is the nearest-neighbor 
intralayer hopping and a = 2.46 A is the graphene lattice 
constant). 

This chiral-invariant Hamiltonian has negative (s = 
— 1) and positive (s = 1) energy eigenstates with eigenen- 
ergies £ Sj fc = S7i|^fc| J and eigenvectors 
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It follows that the interband matrix element of the ve- 
locity operator is 



(s,k\v x \—s,k) = i Jvs\vk\ J 1 sin 



(6) 



Inserting these expressions into the Kubo formula, mul- 
tiplying the result by a factor of g s g v = 4 to account for 
the spin and valley degeneracy of graphene systems, we 
obtain 
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Since 71 is the only energy scale in the Hamiltonian, it is 
clear prior to calculation, that o-r(lo) oc (e 2 /h)(Huj ■ 
The velocity matrix element, joint density-of-states, and 
energy denominator factors combine so that i = for 
every value of J and, importantly, so that ctr(u)) oc J. 
Ideal-model conductivity. — It follows from the chiral 
doublet conductivity Eq. (0), the emergent chiral 
symmetry (l9l l20j of the ideal model, and the chiral- 
ity sum rule Eq. @, that the ideal-model conductivity 
for normal A^-layer graphene satisfies lim^^o o~r(uj) = 



No~ un i. The numerical calculations necessary to evalu- 
ate ctr(w) at intermediate frequencies are also remark- 
ably simple. Because the band energies of the ideal 
model are dependent only on the magnitude of wavevec- 
tor k, the only ^-dependent quantities which appear in 
the wavevector integral for the conductivity (Eq. ([3])) are 
the velocity matrix elements M(<f) = (n,k\hv x \n' ,k). 
In an iV-layer stack the fastest possible angle variation 
in any matrix element varies as exp(±iN(f)). It follows 
that all angle integrals are evaluated exactly by summing 
over 2N + 1 equally spaced orientations. The wavevec- 
tor magnitude integrals are performed by solving for |fc| 
values at which interband energy differences are equal to 
huj. The ideal- model conductivity of any multilayer stack 
can be evaluated very accurately with a relatively small 
numerical effort. 
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FIG. 1: (color online) Conductivity for ideal-model bilayer 
graphene (73 = 0) and for a more realistic model with distant 
neighbor interlayer hopping (73 = 0.3 eV). The inset shows 
ideal-model band structure of bilayer graphene. 



We have evaluated an(uj) curves for many stacking ar- 
rangements; all results confirm our claim that for the 
ideal model o-r(lu) — > Na un i in low and high frequency 
limits. In Figs. HH3] we show representative results for 
some 2, 4, and 10 layer stacks which allow us to discuss 
corrections to the ideal model and trends in the relation- 
ship between stacking and intermediate frequency devia- 
tions from Na un i. 

Discussion. — The optical conductivity of bilayers 
(Fig. [TJ has been studied both theoretically [22I |23| and 
experimentally [13, HH, HE l26l in previous work. The 
electronic structure consists 19] of a J = 2 chiral dou- 
blet at the Fermi energy (the simplest example of emer- 
gent chiral symmetry) and a two-site chain split-off band. 
The low-frequency conductivity originates from transi- 
tions within the J = 2 doublet. Because the correspond- 
ing transition matrix elements are finite at k = 0, the 
strongest infrared (IR) feature occurs at the onset en- 
ergy, Huj = 71, of transitions between the chiral doublet 
and the split-off bands. (71 = 0.3 eV[5] in all our cal- 
culations.) The feature associated with transitions from 
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split-off valence to split-off conduction bands at 271 is 
weaker because in this case the velocity matrix elements 
vanish at the k = onset. Inclusion of remote interlayer 
hopping 73 , or of any other corrections to the ideal model 
has little influence on the optical conductivity. 
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FIG. 2: (color online) Ideal-model band structure and real 
part of the conductivity for all tetralayer graphene stacks, 
ABCA (top), ABCB, ABAC (middle) and ABAB (bottom). 
The insets show stacking diagrams where shaded ovals link 
sublattices a and (5 to the nearest interlayer neighbors. 



The same patterns continue in thicker multilayers. 
The low-frequency conductivity Na un i comes entirely 
from transitions within asymptotically free chiral dou- 
blets, each of which is dispersed across the multilayer. 
The high-frequency conductivity, also Na un i, comes from 
transitions within decoupled single-layers. The crossover 
at intermediate frequencies is punctuated by the onset 
of a series of interband transitions, with the strongest 
features coming from transitions between the low-energy 
chiral doublets and split-off bands. The N — 4 case, il- 
lustrated in Fig. [51 has four distinct stacking sequences, 
two of which are related by inversion symmetry. The 
three inequivalent cases are orthorhombic ABCA stack- 
ing which yields [l^l a J = 4 low-energy chiral doublet 
and three two- site- chain split-off bands, Bernal ABAB 
stacking which yields [l9| two J = 2 chiral doublets and 
four-site-chain split-off bands, and intermediate ABCB 
stacking which yields 19] J — 3 and J — 1 chiral doublets 
and both three and two-site-chain split-off bands. The 
optical conductivity of the orthorhombic ABCA stack 
has a divergent IR feature associated with J = 4 chi- 
ral doublet to two-site chain transitions. The onset of 
this absorption band is below 71 because one of the 



three two-site-chain split-off bands disperses toward the 
Dirac point and has an extremum at finite ka ss 0.1, im- 
plying a divergent joint density of states. The Bernal 
ABAB stack has two jump-discontinuity IR features as- 
sociated with k — transitions between the J = 2 dou- 
blets and E r = ±271 cos(r7r/5) four-site-chain[19] split- 
off bands with r = 1,2. Similarly the ABCB stack has 
strong IR features associated with transitions between 
the chiral doublets and both E = ±71 two-site chain and 
E = ±V27i three-site chain bands. 

Normal graphene stacks can always be organized into 
Bernal and orthorhombic segments. In general stacks 
with more orthorhombic segments have fewer [l9j low- 
energy chiral doublet bands (which must therefore have 
higher chirality because of the chirality sum rule) and 
shorter[19( chain split-off bands. Level repulsion among 
the short-chain bands tends to cause some to disperse to- 
ward the Dirac point and have finite k extrema. Stacks 
with more Bernal segments will have more chiral doublets 
with lower chirality and longer-chain split-off bands. In 
all cases the strongest IR features are associated with 
transitions between chiral doublets and split-off bands. 
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FIG. 3: (color online) Real part of the ideal-model conduc- 
tivity for Bernal (AB), orthorhombic (ABC), and hexagonal 
(AA) 10 layer graphene stacks. 



Since long-chain states are spread over a broader range 
of energies, stacks with more Bernal character tend to 
have IR features that are weaker and spread over a wider 
energy range. This trend is clear in Fig. [3l which com- 
pares the optical conductivities of 10-layer Bernal and 
orthorhombic stacks. In the Bernal case 

[13 

transitions 

from chiral doublets to 10-site-chain split off bands lead 
to jump discontinuities in ctr(u)) at hio = 271 cos(r7r/ll) 
for r = 1,2, 3, 4, 5. In the orthorhombic case, the first IR 
feature appears at a larger frequency but the deviations 
from 10a un i are larger. 

Figure [3] also shows conductivity results for a 10 layer 
hexagonal stack, in which honeycomb layers are placed 
one on top of the other. Stacks of this type do not satisfy 
the chirality sum rule(20j and invariably lead to <7r(ijj) 
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curves which are suppressed at low frequencies. As illus- 
trated in Fig. [31 even N hexagonal stacks do not have a 
low-energy chiral doublet, thus <Jr{uj) vanishes at small 
u) and then increases toward Na un i in steps of 2a un i- 
For odd N the stack has a single J = 1 chiral doublet 
at low-energies so that <jr(u>) starts at a U ni in the low- 
frequency limit and then increases in 2a un i steps toward 
Na un i . More generally, stacks with AA segments always 
have ur(uj) in the low-frequency limit smaller than Na un i 
and approach this limit only at high-frequencies. 

The full electronic structure of graphcnc multilayers is 
usually discussed in terms of an appropriate adaptation 
of the Slonczewski-Weiss-McClure (SWM) [1 S|] param- 
eterization of graphite's bands. A SWM-type model im- 
proves the ideal model by accounting for small differences 
between the energies of 7r-orbitals on inequivalent carbon 
atoms and for various distant neighbor hopping ampli- 
tudes. In the bilayer case the only important additional 
parameter is the distant neighbor interlayer hopping am- 
plitude 73 , and as noted earlier this process has negligible 
influence on ur{ui). These refinements of the electronic 
structure model lead in general to two-dimensional (2D) 
electron and hole Fermi surfaces with Fermi energies that 
are much smaller than the dominant interlayer coupling 
energy 71. The appearance of 2D Fermi surfaces there- 
fore has little influence on the IR conductivity. The pres- 
ence of 2D Fermi surfaces (and in the limit of graphite 
of 3D Fermi surfaces) does imply that or{uj) will in gen- 
eral have a small amplitude Drude peak, not accounted 
for in the present discussion. The Drude peak will take 
a small amount of spectral weight [2!| from the IR inter- 



band transitions. More realistic models also do in general 
break particle-hole symmetry This refinement will cause 
an IR feature associated with a particular split-off va- 
lence band to chiral doublet transition to appear at a 
slightly different frequency than the corresponding tran- 
sition in the ideal model. These caveats notwithstand- 
ing, the origin of the generic <Jr{lo) ~ Na un i behavior 
in graphene multilayers is explained most succinctly by 
the ideal model, and in particular by its emergent chiral 
symmetry. 

The ideal model is also able to capture the stack- 
ing structure implications of measured IR conductiv- 
ity features: the high-frequency decoupled layer Na un i 
limit is approached for hw > 271 w 0.6 eV for nor- 
mal stacks, but only at higher frequencies ~ 471 « 
1.2 eV when AA stacking faults are present. AA 
stacking is also indicated by suppressed conductivity 
at lower frequencies, as discussed previously by Kuz- 
mcnko et al. for the bulk graphite case UJ. In nor- 
mal stacks, more pronounced IR features are an indica- 
tor for orthorhombically stacked sub-units. We conclude 
that the optical conductivity or corresponding transmit- 
tance T(u) = [l + — ctr(u>)] [l8[ spectrum can pro- 
vide a convenient qualitative characterization of multi- 
layer graphene stacks. 
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